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We study thermal and chemical equilibration in 'infinite' hadron matter as well as in 

finite size relativistic nucleus-nucleus collisions using a BUU cascade transport model 

with resonance and string degrees-of-freedom. The 'infinite' hadron matter is simulated 

within a cubic box employing periodic boundary conditions. The various equilibration 

^ ! times depend on baryon density and energy density and are much shorter for particles 

consisting of light quarks then for particles including strangeness. For kaons and antikaons 

O ' the chemical equilibration time is found to be larger than ~ 40 fm/c for all baryon and 

energy densities considered. The inclusion of continuum excitations, i.e. hadron 'strings', 

O ■ leads to a limiting temperature of Tg ~ 150 MeV. 
O 



X 



1. INTRODUCTION 

O I Nucleus-nucleus collisions at relativistic and ultrarelativistic energies are studied exper- 
^ ' imentally and theoretically to obtain information about the properties of hadrons at high 
density and/or temperature as well as about the phase transition to a new state of matter, 
the quark-gluon plasma (QGP). In the latter deconfined partons are the essential degrees 
of freedom that resolve the underlying structure of hadrons |]1| . Whereas the early 'big- 
I bang' of the universe most likely evolved through steps of kinetic and chemical equilibrium, 
the laboratory 'tiny bangs' proceed through phase-space configurations that initially are 
far from an equilibrium phase and then evolve by fast expansion. These 'specific initial 
conditions' - on the theoretical side - have lead to a rapid development of nonequilibrium 
quantum field theory and nonequilibribrium kinetic theory Presently, semiclassical 

transport models are widely used as approximate solutions to these theories and practi- 
cally are an essential ingredient in the experimental data analysis. For recent reviews we 
refer the reader to Refs. [^H^. 

On the other hand, many observables from strongly interacting systems are dominated 
by many-body phase space such that spectra and abundances look 'thermal'. It is thus 
tempting to characterize the experimental observables by global thermodynamical quan- 
tities like 'temperature', chemical potentials or entropy [[^-[TT[|. We note, that even the 
use of macroscopic models like hydrodynamics [IT^IO] employs as basic assumption the 



concept of local thermal and chemical equilibrium. The crucial question, however, how 
and on what timescales a global thermodynamic equilibrium can be achieved, is presently 
a matter of debate. Thus nonequilibrium approaches have been used in the past to ad- 
dress the problem of timescales associated to global or local equilibration fl^-^. In 
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view of the increasing 'popularity' of thermodynamic analyses a thorough microscopic 
reanalysis of this questions appears necessary especially for nucleus-nucleus collisions at 
ultrarelativistic energies that aim at 'detecting' a phase transition to the QGP. 

In this contribution we discuss equilibration phenomena in 'infinite' hadronic matter 
using a microscopic transport model that contains both hadron resonance and string 
degrees-of-freedom. With this investigation we want to provide insight into the reaction 
dynamics by the use of cascade-like models and also point out some of their limitations. 
The 'infinite' hadronic matter is modelled by initializing the system solely by nucleonic 
degrees of freedom through a fixed baryon density and energy density, while confining 
it to a cubic box and imposing periodic boundary conditions during the propagation in 
time. 



2. EQUILIBRATION AND LIMITING TEMPERATURE 

To investigate the equilibration phenomena addressed above we perform microscopic 



calculations using the Boltzmann-Uehling-Uhlenbeck (BUU) model of Refs. p3|,^. This 
model is based on the resonance concept of nucleon-nucleon and meson-nucleon interac- 
tions at low invariant energy ^/s ||2^, adopting all resonance parameters from the Manley 



analysis [^. The high energy collisions - above a/s = 2.6 GeV for baryon-baryon colli- 
sions and a/s = 2.2 GeV for meson-baryon collisions - are described by the LUND string 
fragmentation model FRITIOF [Q. This aspect is similar to that used in the HSD ap- 
proach |]6|j27H29| and the UrQMD code 0. For a detailed description of the underlying 



model at low energy we refer the reader to Ref. P^ |. 

2.1. A box with periodic boundary conditions 

In order to study 'infinite' hadronic matter problems we confine the particles in a 
cubic box with periodic boundary conditions for their propagation similar to a recent 



box calculation within the UrQMD model |T^. We specify the initial conditions, i.e 



baryon density p, strange particle density ps and energy density e as follows: first the 
initial system is fixed to Np = 80 protons and Nn = 80 neutrons, which are randomly 
distributed in a cubic box of volume V. The 3-momenta pi of the nucleons in a first step 
are randomly distributed inside a Fermi-sphere of radius pp = 0.26 GeV/c (at po) and 
in a second step boosted by ±/5cm by a proper Lorentz transformation. Thus the initial 
baryon density p is fixed as p = A/V, A = Np + Nn- The strange particle density is set to 
zero as in related heavy-ion experiments while the energy density is defined as e = E/V, 
where E is the total energy of all nucleons 

A 



E = Y: + m%. (1) 

i 

The boost velocity (3cm is related to the initial energy density e (excluding Fermi motion) 
as 



P^ = ^jl-'-^ (2) 



using e = •jcmP^N with ■jcm = l/y'l — Pcm- Recall that pofriN ~ 0.15 GeV/fm^ so that 
an energy density e ~ 1.5 GeV/fm^ at density po corresponds to — 10, i.e. the SPS 
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energy Tiab — 185 A-GeV. We thus start with a 'true' nonequihbrium situation in order to 
mimique the initial stage in a relativistic heavy-ion colhsion. The initial phase represents 
two interpenetrating, (ideally) infinitely extended fiuids of cold nuclear matter. 

We now propagate all particles in the box in the cascade mode (without mean-field 
potentials) using periodic boundary conditions, i.e. particles moving out of the box are 
reinserted at the opposite side with the same momentum. The phase-space distribution 
of particles then can change due to elastic collisions, resonance and string production and 
their decays to mesons and baryons again. We recall that we include all baryon resonances 
up to an invariant mass of 2 GeV and meson resonances up to the 0-meson. According 
to the initial conditions for e and p the factor determines if strings are excited in the 
very first collisions. This is the case for > 1.4 where the early equlibration stages are 
dominated by string formation and decay. 

2.2. Chemical equilibration 

Figure |I| shows the time evolution of the various particle abundances (nucleons A^, 
A, A, TT and mesons) for density p = po (left panel) at energy density e = 0.22 
GeV/fm^ and for p = 3po (right panel) at e = 0.66 GeV/fm^. These initial conditions 
correspond to bombarding energy Tiab per nucleon of roughly 2 A-GeV. In Fig. |1| we count 
all particles which are 'hadronized', i.e. produced by string decay after a formation time 
of Tp = 0.8 fm/c in their rest frame. 
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Figure 1. Time evolution of the various particle abundances (nucleons A^, A, A, vr and 
mesons) for density p = po (left panel) at energy density e = 0.22 GeV/fm^ and for 
p = 3po (right panel) at e = 0.66 GeV/fm^. 

After several fm/c the number of nucleons decreases due to inelastic collisions that 
produce either baryon resonances or additional mesons. The number of A-resonances 
grows up to a maximum in a few fm/c, since a lot of A's are produced in the first A^A^ 
collisions; their number subsequently decreases with time due to their decay and excitation 
of further resonances or due to reabsorption. The numbers of vr's and t^'s increase very 
fast and reach the equilibrium value within a few fm/c whereas the strange particles 
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{K^, K~,A) require a much longer time for equilibration. 

For the higher energies the initial particle production proceeds via the formation and 
decay of string excitations. This leads in particular to a very early onset of strange 
particles (mainly kaons and hyperons) within the first fm/c either due to the initial strings 
or due to secondary or ternary baryon-baryon, meson-baryon and meson-meson induced 
string-like interactions. In Ref. it was shown that these early secondary and ternary 
reactions can contribute up to about 50 % of the total strange particles obtained in a 
Pb + Pb reaction at CERN SPS energies and thus explain the factor of 2 in the observed 
relative strangeness enhancement compared to p+p reactions. This, however, does not 
imply that chemical equilibrium for the dominant strange particles has been achieved in 
this reaction, as our analysis clearly shows. In the later stages, when the system has 
become, more or less, isotropic in momentum space, strange particles can only be further 
produced by low energy hadronic reactions, which, however, have a considerable threshold 
and are thus strongly suppressed. This explains the long chemical equilibration times for 
the strange particles first demonstrated by Koch, Miiller and Rafelski |]T^ . 



In order to define an overall chemical equilibration time we perform a fit to the particle 
abundances N{t) for pions and kaons as 

iV(t) =iVe,(l-exp(-t/re,)), (3) 

where N^q is the equilibrium limit. The equilibration time Teg thus corresponds to the 
time t when ~ 63% of Neq is achieved. 

Figure |^ shows the equilibration time r^g versus energy density for tt and mesons 
at different baryon densities of l/3po5PO)3po and 6po- We find that the equilibration 
time for pions scales as rj^^ ~ 1/p or T^^ ~ p, thus we present the curve only for baryon 
density po. Whereas rj^ slowly grows with energy-density, falls steeply with e. This 
marked difference is due to the fact that, on one hand, the kaon production rate increases 
dramatically with ^/s whereas that of the pions, on the other hand, is more fiat. With 
increasing energy thus more strange particles are produced through strings especially from 
the primary collisions with high ^/s and the chemical equilibration is achieved faster. 

In Fig. 1^ we have considered an 'ideal' situation, i.e. hadron matter at fixed energy 
and baryon density. In realistic heavy-ion collisions the system goes through the different 
stages due to interactions and expansion. However, as follows from Fig. ^, the equili- 
bration time for strangeness is larger than 40 fm/c for all energy and baryon densities. 
Thus in realistic nucleus-nucleus collisions the chemical equilibration of strange particles 
requires also a time above 40 fm/c which is considerably larger than the actual reaction 
time of a few 10 fm/c or less pO| . 

The particle abundances used to extract T^g in Fig. ^ have been calculated without any 
in- medium potentials. In fact, the introduction of attractive potentials (especially for K~) 
will lower the hadronic thresholds and thus increase the scattering rate between strange 
and nonstrange hadrons, whereas the feels some repulsive potential and the trend 
goes in the opposite way. According to our calculations such in-medium modifications (in 
line with Ref. |^) give a correction to the equilibration times by atmost 10 % and 
shortens the K~ equilibration times up to 20 % at density po- 
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Figure 2. Equilibration time Teq versus energy density e for vr and mesons at different 
baryon densities l/Spo^Po^Spo and 6po- 



2.3. Thermal equilibration and limiting temperature 

In this subsection we investigate the approach to thermal equilibration. For the equi- 
librated system we can extract a temperature T by fitting the particle spectra with the 
Boltzmann distribution 



exp{-Ei/T), 



(4) 



where Ei = yPi+m^ is the energy of particle i. We note that at the temperatures of 
interest here, the Bose and Fermi distributions are practically identical to a Boltzmann 
distribution. We find that in equilibrium the spectra of all particles can be characterized 
by one single temperature T |^ . 

In the left panel of Fig. ^ we display the temperature T versus energy density e for 
different baryon densities p: l/3po (open down triangles), po (full squares), 3po (full dots), 
6po (full up triangles). In order to compare calculations for different baryon densities we 
have subtracted the baryon energy density at rest, i.e. ~ m^rp (except for Fermi motion). 
As seen from Fig. ^ the temperature grows with energy density up to a limiting value 



reminiscent of a 'Hagedorn' temperature 0]. From our detailed investigations we obtain 
for the limiting temperature T, ~ 150 ± 5 MeV which practically does not depend on 
baryon density. Such a singular behavior of s{T) for T c::^ Tg has also been found in the 
box calculations in Ref. 



for p = pq. Our limiting temperature is slightly higher than 
that in Ref. [|l8l {Tg = 130±10 MeV) due to the different number of degrees of freedom; the 
model [|I^] contains more resonances and uses a different threshold for string excitations. 
Thus, there is some phenomenological sensitivity to the hadronic zoo of particles and 
string thresholds employed in the model. 

In order to investigate the equilibrium behavior of hadron matter we also compare our 
transport (box) calculations with a simple Statistical Model (SM) for an Ideal Hadron Gas 
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Figure 3. Left panel: equilibrium temperature T versus the energy density e — rriNp for 
different baryon densities p: l/3po (open down triangles), po (foU squares), 3po (full dots), 
6po (full up triangles). Right panel: equilibrium temperature T versus the energy density 
for baryon density p = po- The full dots correspond to the statistical model (SM) without 
strings, the full squares show our box calculations including string degrees of freedom, 
while the solid line shows the result from the extended SM including a Hagedorn mass 
spectrum for strings. 



where the system is described by a grand canonical ensemble of non-interacting fermions 
and bosons in equilibrium at temperature T. All baryon and meson species considered in 
the transport model also have been included in the statistical model ||30(| . 



Within the SM we find that the temperature increases continuously with energy density 
since the continuum excitations, i.e. the string degrees of freedom, are not included (full 
dots in the right panel of Fig. ^), whereas the box calculation with strings gives the 
limiting temperature (full squares in Fig. |^). Both curves in Fig. ^ have been calculated 
for density po. 

To reproduce qualitatively our box result within the SM we have to include continuum 
excitations in the statistical model, i.e. a Hagedorn mass spectrum for strings |^ (for 
details see [0). For the 'Hagedorn' temperature we use the temperature Tg as ob- 
tained from the box calculations, i.e. Th = ~ 150 MeV. As seen in right panel of Fig. ^ 
we achieve agreement of the extended SM and our box calculations. However, we point 
out that the limiting temperature Tg from our string model involves somewhat different 
physics assumptions than the Hagedorn model at temperature Th- Tg should not really 
be identified with the 'Hagedorn' temperature Th, though close similarities exist. In the 
Hagedorn picture and for temperatures close to Th the abundance of 'normal' hadrons 
or known resonances stays constant with increasing energy density whereas the number 
and energy density of the (hypothetical) bootstrap excitations diverges for T Th- The 
Hagedorn model thus assumes 'particles' of mass m — > oo to be populated for T — > Th, 
that dynamically can be formed in collisions of high mass hadrons for t — oo. In contrast, 
our string model does not include energetic string-string interactions that might produce 
more massive strings. 
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3. SUMMARY 

In this contribution we have performed a systematic study of equihbration phenom- 
ena and equihbrium properties of 'infinite' hadronic matter as well as of relativistic 
nucleus-nucleus collisions using a BUU transport model that contains resonance and string 
degrees-of-freedom. The 'infinite' hadron matter is modelled by initializing the system at 
fixed baryon density, strange density and energy density by confining it in a cubic box 
with periodic boundary conditions ||30|| . 

We have shown that the equilibration times Tgg for different particles depend on baryon 
density and energy density. The time Tf.q for non-strange particles is much shorter than 
for particles including strangeness; for kaons and antikaons the equilibration time is found 
to be larger than ~ 40 fm/c for all baryon and energy densities considered. The overall 
abundance of the dominant strange particles (kaons and A's) being produced and obtained 
within the BUU cascade model for heavy-ion collisions can therefore not be described by 
assuming a perfect chemical equilibrium as strangeness is typically still undersaturated to 
a quite large extent. We mention that transport model calculations like ours can describe 
the yield and spectra of the produced nonstrange hadrons as well as ^ K~ ,A yields quite 
well at SPS energies [|6|,p9|. On the other hand, at AGS energies the measured K^/tt~ 



ratio in central Au + Au collisions is underestimated by about 30% [^. However, we 
have to point out that the more exotic strange particles (like the measured antihyperon 
yields of Ref. [^) can by far not be explained within such standard hadronic multiple 
channel reactions. These hadronic data possibly point towards new physics. 

We have, furthermore, shown that thermal equilibrium is established quickly, within 
about 5 fm/c at SIS energies and samewhat larger times at high energies. The inclusion 
of continuum excitations, i.e. hadron 'strings', leads to a limiting temperature of Tg ~ 
150 MeV in our transport approach which practically does not depend on the baryon 
density and energy. We have compared our results with the statistical model (SM), which 
contains the same degrees of freedom and the same spectral functions of particles as our 
transport model. We found that the limiting temperature behaviour can be reproduced 
in the statistical model only after including continuum excitations of the Hagedorn type, 
otherwise the fireball temperature extracted from the particle abundances and spectra is 
overestimated substantially. 

Close to the critical temperature T^, the hadronic energy densities can increase to a 
couple of GeV/fm^. From lattice QCD calculations one expects that a phase transition to 
a potentially deconfined QGP state should occur. Referring to the limiting temperature 
Tg ^ 150 MeV obtained, a QGP should be revealed and clearly distinguished from a 
hadronic state of matter if one can unambiguously prove the existence of an equilibrated 
and thermal phase of strongly interacting matter with temperatures exceeding, e.g., 200 
MeV. The best candidates are electromagnetic probes, either direct photons or dileptons. 
On the other hand these are also 'contaminated' by hadronic background and/or pree- 
quilibrium physics. So far no thermal electromagnetic source with temperatures larger or 
equal than 200 MeV has been clearly identified. This might happen at RHIC energies in 
central Au + Au collisions which are expected to be studied soon. 
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